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Abstract 

Given two polygonal curves in the plane, there are many ways to define a notion of similarity between them. 
One measure that is extremely popular is the Frechet distance. Since it has been proposed by Alt and Godau 
in 1992, many variants and extensions have been studied. Nonetheless, even more than 20 years later, the 
original 0(n 2 log n) algorithm by Alt and Godau for computing the Frechet distance remains the state of 
the art (here n denotes the number of vertices on each curve). This has led Helmut Alt to conjecture that 
the associated decision problem is 3SUM-hard. 

In recent work, Agarwal et al. show how to break the quadratic barrier for the discrete version of 
the Frechet distance, where one considers sequences of points instead of polygonal curves. Building on their 
work, we give a randomized algorithm to compute the Frechet distance between two polygonal curves in time 
0(n 2 v / logn(log log n) 3 / 2 ) on a pointer machine and in time 0(n 2 (loglogn) 2 ) on a word RAM. Furthermore, 
we show that there exists an algebraic decision tree for the decision problem of depth 0(n 2 ~ 6 ), for some 
e > 0. This provides evidence that the decision problem may not be 3SUM-hard after all and reveals an 
intriguing new aspect of this well-studied problem. 



1 Introduction 

Shape matching is a fundamental problem in computational geometry, computer vision, and image 
processing. A simple version can be stated as follows: given a database T> of shapes (or images) 
and a query shape S, find the shape in T> that most resembles S. However, before we can solve 
this problem, we first need to address a much more fundamental issue: what does it mean for two 
shapes to be similar? In the mathematical literature, there are many different notions of distance 
between two sets, a prominent example being the Hausdorff distance. For two planar sets A and 
B, the (undirected) Hausdorff distance is denned as 

dj{{A, B) := maxjsup inf \\a — b\\, sup inf \\a — b\\}, 

where || • || denotes the Euclidean distance. It has the advantage of being simple to describe and easy 
to compute for discrete sets. In the context of shape matching, however, the Hausdorff distance 
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Fig. 1: Two polygonal curves with low HausdorfT distance but high Frechet distance. 



often turns out to be unsatisfactory. There are well known examples where the distance fails to 
capture the similarity of shapes as perceived by human observers (2). The problem is that it does 
not take the continuity of the shapes into account. Figure [T] shows an example. 

In order to address this issue, Alt and Godau introduced the Frechet distance into the com- 
putational geometry literature [6,33,. They argue that the Frechet distance is better suited as a 
similarity measure, and they describe an 0(n 2 logn) time algorithm to compute it on a real RAM 
or pointer machine^] Since Alt and Godau's seminal paper, there has been a wealth of research 
in various directions, such as extensions to higher dimensions (5 17 19 
tion algorithms 
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38 - |40] . All known approximation algorithms make further assumptions on the 
curves, and no o(n 2 log n)-time approximation algorithm is known for arbitrary polygonal curves. 
The Frechet distance and its variants, such as dynamic time- warping |10| , have found various ap- 
plications, with recent work particularly focusing on geographic applications such as map-matching 
tracking data [IT, 43 and moving objects analysis (13 14,35). 

Despite the large amount of published research, the original algorithm by Alt and Godau has not 
been improved, and the quadratic barrier on the running time of the associated decision problem 
remains unbroken. If we cannot improve on a quadratic bound for a geometric problem despite 
many efforts, a possible culprit may be the underlying 3SUM-hardness (32). This situation induced 
Helmut Alt to make the following conjecture]^] 

Conjecture 1.1 (Alt's Conjecture). Let P and Q be two polygonal curves in the plane. Then 
it is 3S\JM-hard to decide whether the Frechet distance between P and Q is at most 1. 

Here, 1 can be considered as an arbitrary constant, which can be changed to any other bound 
by scaling the curves. So far, we know only that this problem takes f2(nlogn) steps in the algebraic 
computation tree model (15) . 

Recently, Agarwal et al. [l] showed how to achieve a subquadratic running time for the discrete 
version of the Frechet distance. Their approach relies on reusing small parts of the solution. 
We follow a similar approach based on the so-called Four-Russian-trick which precomputes small 
recurring parts of the solution and uses table-lookup to speed up the whole computation^] The 
result by Agarwal et al. is stated in the word RAM model of computation. They ask whether their 
result can be generalized to the case of the original (continuous) Frechet distance. 

Our contribution. We address the question by Agarwal et al. and show how to extend their 
approach to the Frechet distance between two polygonal curves. Our algorithm requires total 
expected time 0{n 2 \/log n(log log n) 3 / 2 ). This is the first algorithm to achieve a running time of 



1 For a brief overview of the different computational models in this paper, refer to Appendix |X] 

2 Personal communication 2012, see also [5]. 

3 It is well known that the four Russians are not actually Russian, so we refer to them as four Soviets in the title. 
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o(n 2 log n) and constitutes the first improvement for the general case since the original paper by 
Alt and Godau [6]. We emphasize that our algorithm runs on a real RAM/pointer machine and 
does not require any bit manipulation tricks to achieve the speedup. If we relax the model to 
allow constant time table-lookups, the running time can be improved to be almost quadratic, up to 
O (log log n) factors. As in Agarwal et at, our results are achieved by first giving a faster algorithm 
for the decision version, and then performing an appropriate search over the critical values to solve 
the optimization problem. 

In addition, we show that non-uniformly, the Frechet distance can be computed in subquadratic 
time. More precisely, we prove that the decision version of the problem can be solved by an algebraic 
decision tree (9] of depth 0(n 2_E ), for some fixed e > 0. This makes it unlikely that the Frechet 
distance is 3SUM-hard, since it is conjectured that no such decision tree exists for 3SUM [3]. 
However, it is not clear how to implement this decision tree in subquadratic time, which hints at 
a discrepancy between the decision tree and the uniform complexity of the Frechet problem. This 



puts it into the illustrious company of such notorious problems as Sorting X + Y 31 , Min-Plus- 



Convolution [12] , or finding the Delaunay triangulation for a point set that has been sorted in 
two orthogonal directions [20] . We find that this aspect of the Frechet distance is highly intriguing 
and deserves further study. 



2 Preliminaries and Basic Definitions 



Let P and Q be two polygonal curves in the plane, defined by their vertices PojPl> • • • iPn an d 
qo, qx,..., q n ^\ Depending on the context, we interpret P and Q either as sequences of n edges, or 
as continuous functions P, Q : [0, n] — > M 2 . In the latter case, we have P(i + A) = (1 — \)pi + \p%+i 
for i = 0, ... ,n — 1 and A G [0, 1], and similarly for Q. Let \£ be the set of all continuous and 
nondecreasing functions a: [0, n] — > [0, n] with <x(0) = and cr(n) = n. The Frechet distance 
between P and Q is defined as 

d F (P,Q) := inf max \\P(x) - Q(a(x))\\, 
o-e* xe[o,n] 

where || • || denotes the Euclidean distance. 

The free-space diagram. The classic approach to compute Q) uses the free-space diagram 

FSD(P,Q). It is defined as 

FSD(P,Q) := {(x,y) G [0,n] 2 | \\P(x) - Q(y)\\ < 1}. 

In other words, FSD(P, Q) is the subset of the joint parameter space for P and Q where the 
corresponding points on the curves have distance at most 1, see Figure [2] 

The structure of FSD(P, Q) is easy to describe. Let R := [0, n] 2 be the ground set. We subdivide 
R into n 2 cells C(i,j) = + x [j,j + 1], for i,j = 0, ...,n— 1. The cell C(i,j) corresponds 
to the edge pair a+i and fj+i, where e^+i is the (i + l) th edge of P and fj+i is the (j + l) th 
edge of Q. Then the set F(i,j) := FSD(P, Q) n C(i,j) represents all pairs of points on e^+i x 
with distance at most 1. Elementary geometry shows that F(i,j) is the intersection of C(i,j) with 
an ellipse pi. In particular, the set F(i,j) is convex, and the intersection of FSD(P, Q) with the 

4 For simplicity, we assume that both curves have the same number of vertices. This constitutes the most interesting 
case from a theoretical point of view, and it is straightforward to extend our results to the case when the number of 
vertices on the two curves differs. 
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Fig. 2: Two polygonal curves P and Q and their associated free-space diagram. For example, the 
white area in C(2, 1), denoted P(2, 1), corresponds to all points on the third edge of P and 
the second edge of Q that have distance at most 1. The region reach(P, Q) is shown in blue. 
Since (5, 5) is reachable from (0, 0), we have cIf(P, Q) < 1. 

boundary of C(i,j) consists of four (possibly empty) intervals, one on each side of dC(i,j). We 
call these intervals the doors of C(i,j) in FSD(P, Q). A door is said to be closed if the interval is 
empty, and open otherwise. 

A path 7r in FSD(P, Q) is bimonotone if it is both x- and y-monotone, i.e., every vertical and 
every horizontal line intersects it in at most one connected component. Alt and Godau observed 
that it suffices to determine whether there exists a bimonotone path from (0, 0) to (n, n) inside 
FSD(P, Q). More precisely, define 

reach(P, Q) := {p G FSD(P, Q) \ p is reachable from (0, 0) in FSD(P, Q) on a bimonotone path}. 

Then dp(P,Q) < 1 if and only if (n,n) £ reach(P, Q). It is not necessary to compute all of 
reach(P, Q): since FSD(P, Q) is convex inside each cell, we actually need just the intersections 
reach ( P, Q) n dC(i,j). The sets defined by reach ( P, Q) n dC(i,j) are subintervals of the doors of 
the free-space diagram, and they are defined by endpoints of doors in the free-space diagram in 
the same row or column. We call the intersection of a door with reach(P, Q) a reach-door. The 
intersections can be found in 0(n 2 ) time through a simple traversal of the cells jjj]. In the next 
sections, we show how to obtain the crucial information, namely whether (n,n) £ reach(P, Q), in 
o(n 2 ) instead. 

Basic approach and intuition. In our algorithm for the decision problem, we basically want 
to compute reach(P, Q). But instead of propagating the reachability information cell by cell, we 
always group r by r cells (with 1 <C r <C n) into an elementary box of cells. When processing 
a box, we can assume that we know which parts of the left and the bottom boundary of the box 
are reachable. That is, we know the reach-doors on the bottom and left boundary, and we need to 
compute the reach-doors on the top and right boundary of the elementary box. 

The reach-doors on the top and right are determined by the combinatorial structure of the box. 
More specifically, if we know for every row and column the order of the door endpoints (including 
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the reach-doors on the left and bottom boundary), we know which door boundaries determine the 
reach-doors on the top and right boundary. We will call the sequence of these orders, the (full) 
signature of the box. 

The total number of possible signatures is bounded by an expression in r. Thus, if we pick r 
sufficiently small compared to n, we can pre-compute for all possible signatures the reach-doors on 
the top and right boundary, and build a data structure to query these quickly (Section [3]). Since 
the reach-doors on the bottom and left boundary are required to make the signature, we initially 
have only partial signatures. In Section |4j we describe how to compute these efficiently. The partial 
signatures are then used to preprocess the data structure such that we can quickly find the full 
signature once we know the reach-doors of an elementary box. After building and preprocessing 
the data structure, it is possible to determine oIf(P,Q) < 1 efficiently by traversing the free-space 
diagram elementary box by elementary box, as explained in Section [5] 

3 Building a lookup table 

3.1 Preprocessing an elementary box 

Before even considering the input, our algorithm builds a lookup table. As mentioned above, the 
purpose of this table is to speed up the computation of small parts of the free-space diagram. 

Let r £ N be a parameter^] The elementary box is a subdivision of [0,r] 2 into r columns and 
rows, thus t 2 cellsj^] For i,j = 0, ... ,r — 1, we denote the cell + 1] x [j,j + 1] with D(i,j). 
We denote the left side of the boundary dD(i,j) by l(i,j) and the bottom side by b(i,j). Note 
that l(i,j) coincides with the right side of dD(i — 1, j) and that b(i,j) coincides with the top of 
dD(i,j — 1). Thus, we write 1(t + for the right side of D(r,j) and b(i,r + 1) for the top side 
of D(i,r). Figure [3] shows the elementary box. 

5 A preview for the impatient reader: we will later set r = 0(\/logn/ log log n). 

6 For now, the elementary box is a combinatorial concept. In the next section, we overlay these boxes on the 
free-space diagram to obtain "concrete" elementary boxes. 

Ti 1 I 1 1 



i(2,l) 



b(2, 2) 

D(2,l] 

b(2, 1) 



1(3.1) 



Fig. 3: The elementary box. The cell D(2, 1) is shown in white. The left vertical boundaries 1(2, 1) 
and 1(3, 1) are red, the horizontal boundaries 6(2, 1) and 6(2, 2) are orange. 
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Fig. 4: The door-order of a row encodes the combinatorial structure of the doors. The door-order 
for the row in the figure is siS3S4tst3toS2t4SoS5tit2. Note that so an d to represent the 
reach-door, which is empty in this case. These are omitted in the partial door-order. 



The door-order crj for a row j is a permutation of {so, to, ... , s T , t T }, thus having 2r+2 elements. 
For i = 1, . . . ,t, the element Si represents the lower endpoint of the door on l(i, j), and U represents 
the upper endpoint. The elements sq and to are an exception: they describe the reach-door on the 
left boundary 1(0, j) (that is, its intersection with reach(P, Q)). The door-order <rj represents the 
combinatorial order of these endpoints, as projected onto a vertical line (see Figure [4b. 

Note that some door-orders may encode the same combinatorial structure. In particular when 
door i is closed, the exact position of Sj and tj in a door-order is irrelevant, up to tj being before 
Si. For a closed door i (i > 0), we assign Sj to the upper endpoint of and tj to the lower 

endpoint. The values of so and to are defined by the reach-door and their relative order is thus a 
result of computation. We break ties between Sj and by placing Sj before tj/, and any other ties 
are resolved by index. A door-order af is defined analogously for a column i. We write x <\ y if x 
comes before y in af, and x < r j y if x comes before y in a r -. A partial door-order is a door-order in 
which so and to are omitted (i.e. the intersection of reach(P, Q) with the door is still unknown). 

We define the (full) signature of the elementary box as the aggregation of the door-orders of 
its rows and columns. Hence, a signature S = (<rf , . . . , a\, . . . , a r T ) consists of 2r door-orders: 
one door-order a\ for each column % and one door-order a r ^ for each row j of the elementary box. 
Similarly, a partial signature is the aggregation of partial door-orders. 

For a given signature, we define the combinatorial reachability structure of the elementary box 
as follows. For each column i and for each row j, the combinatorial reachability structure indicates 
which door boundaries in the respective column or row define the reach-door of b(i, r) or l(r,j). 

Lemma 3.1. Let £ be a signature for the elementary box. Then we can determine the combinatorial 
reachability structure o/E in total time 0(r 2 ). 

Proof. We use dynamic programming, very similar to the algorithm by Alt and Godau [6] . For each 
vertical edge l(i,j) we define a variable l(i,j), and for each horizontal edge b(i,j) we define a variable 
b{i,j). The l(i,j) are pairs of the form (s u ,t v ), representing the reach-door reach ( P, Q) n l(i,j). 
If this reach-door is closed, then t v < r j s u holds. If the reach-door is open, then it is bounded by 
the lower endpoint of the door on l(u,j) and by the upper endpoint of the door on l(v,j). (Note 
that in this case we have v = i.) Once again sq and to are special and represent the reach-door on 
1(0, j). The variables b(i,j) are defined analogously. 

Now we can compute l(i,j) and b(i,j) recursively as follows: first, we set 

= b(i,0) = (s , t ), for i,j = 0,...,t- 1. 

Next, we describe how to find l(i,j) given l(i — and b(i — 1, j), see Figure^! 



6 




Fig. 5: The three cases for the recursive definition of If the lower boundary is reachable, 

we can reach the whole right door (left). If neither the lower nor the left boundary is 
reachable, the right door is not reachable either (middle). Otherwise, the lower boundary 
is the maximum of l{i — and the lower boundary of the right door (right). 

Case 1: Suppose b(i — is open. This means that b(i — l,j) intersects reach(P, Q), so 
reach(P, Q) n is limited only by the door on and we can set := (sj, U). 

Case 2: If both b(i — l,j) and l(i — 1, j) are closed, it is impossible to reach and thus we set 

T(i,j):=T(i-l,j). 

Case 3: If b(i — is closed and l(i — is open, we may be able to reach via l(i — 

Let s u be the lower endpoint of 1, j). We need to pass above s u and S{ and below ti, and 

therefore set := (max(s M , S{), ti), where the maximum is taken according to the order < r -. 

The recursion for the variable b(i,j) is defined similarly. Now it is easy to see that we can 
implement the recursion in time 0(t 2 ) for any given signature, for example by traversing the 
elementary box column by column, while processing each column from bottom to top. □ 

Clearly, there are at most ((2r + 2)!) 2t = t°( t2 ) distinct signatures for the elementary box. 
We will choose r = A-^/logn/log logn for a sufficiently small constant A > 0, so that this number 
becomes o(n). Thus, during the preprocessing stage we have time to enumerate all possible signa- 
tures and determine the corresponding combinatorial reachability structure inside the elementary 
box. This information is then stored in an appropriate data structure. 

3.2 Building the data structure 

Before we can describe this data structure, we must first explain how the door-orders are repre- 
sented. This depends on the computational model. Note that by our choice of r, there are o{n) 
distinct door-orders. On the word RAM, we can represent each door-order and partial door-order 
by an integer between 1 and (2r)!. This certainly fits into a word of logn bits. On the pointer ma- 
chine, we create a record for each door-order and each partial door-order. Then, we can represent 
an order by a pointer to the corresponding record. 

We now explain how we organize the data structure. It consists of two stages, as schematized 
in Figure [6| In the first stage (Figure [6] a-b), we assume we know the partial door-order for each 
row and for each column of the elementary bosQ and we wish to determine the partial signature. 
In the second stage (Figure [6] c-d), we have obtained the reach-doors for the left and bottom sides 
of the elementary box, and we are looking for the full signature. The details of our method depend 

7 In the next section, we describe how to determine the partial door-orders efficiently. 
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a) partial door-orders 

b) partial signature 

c) reach-doors 

d) full signature 



combinatorial reachability 

Fig. 6: Overview of the data structure, a-b) Using partial door-orders, we find the partial signature 
of an elementary box. c-d) Once we know the reach-doors on the bottom and left boundary 
of an elementary box, we determine the full signature. This results in the combinatorial 
reachability structure of that box. 

on the computational model. One way uses table lookup and requires the word RAM, the other 
way works on the pointer machine, but is a bit more involved. 

Word RAM. We organize the lookup table as a large tree T. In the first stage, each level of 
T corresponds to a row or column of the elementary box. Thus, there are 2r levels. Each node 
has (2t)! children, representing the possible partial door-orders for the next row or column. Since 
we represent door-orders by positive integers, each node of T may store an array for its children, 
and we can choose the appropriate child for a given partial door-order in constant time. Thus, to 
determine the partial signature for an elementary box requires 0(t) steps on a word RAM. 

For the second stage, we again use a tree structure. Now, however, the tree has O(r) layers, 
each with O(logr) levels. Again, each layer corresponds to a row or column of the elementary 
box. The levels inside each layer then implement a balanced binary search tree that allows us to 
locate the endpoints of the reach-door within the partial signature. Since there are 2r endpoints, 
this requires O(logr) levels. Thus, it takes 0(r log r) time to find the full signature of a given 
elementary box. 

Pointer machine. Now we are not allowed to store a lookup table on every level of the tree T, 
and there is no way to quickly find the appropriate child for a given door-order. Instead, we must 
rely on batch processing to achieve a reasonable running time. 

Thus, suppose that during the first stage we want to find the partial signatures for a set B of m 
elementary boxes, where again for each box in B we know the partial door-order for each row and 
each column. Recall that we represent the door-order by a pointer to the corresponding record. 
With each such record, we store a queue of elementary boxes that is empty initially. 

We simultaneously propagate the boxes in B through T, proceeding level by level. In the first 
level, all of B is assigned to the root of T. Then, we go through the nodes of one level of T, from 
left to right. Let v be the current node of T. We consider each elementary box b assigned to 
v. We determine the next partial door-order for b, and we append b to the queue for this partial 
door-order — the queue is addressed through the corresponding record, so all elementary boxes with 
the same next partial door-order end up in the same queue. Next, we go through the nodes of the 
next level, again from left to right. Let v' be the current node. The node v' corresponds to a next 
partial door-order a that extends the known signature of its parents. We consider the queue stored 



root 
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at the record for a. By construction, the elementary boxes that should be assigned to v' appear 
consecutively at the beginning of this queue. We remove these boxes from the queue and assign 
them to v'. After this, all the queues are empty, and we can continue by propagating the boxes to 
the next level. During this procedure, we traverse each node of T a constant number of times, and 
in each level of the T we consider all the boxes in B. Since T has o(n) nodes, the total running 
time is 0{n + tot). 

For the second stage, the data structure works just as in the word RAM case, because no table 
lookup is necessary. Again, we need O(rlogr) steps to process one box. 

After the second stage we obtain the combinatorial reachability structure of the box in constant 
time since we precomputed this information for each box. Thus, we have shown the following 
lemma, independently of the computational model. 

Lemma 3.2. For r = Ay^logn/log logn with a sufficiently small constant A > 0, we can construct 
in o(n) time a data structure of size o{n) such that 

• given a set of to elementary boxes where the partial door-orders are known, we can find the 
partial signature of each box in total time 0{n + mr); 

• given the partial signature and the reach- doors on the bottom and left boundary of an elemen- 
tary box, we can find the full signature in 0(r log r) time; 

• given the full signature of an elementary box, we can find the combinatorial reachability struc- 
ture of the box in constant time. 

□ 

4 Preprocessing a given input 

Next, we perform a second preprocessing phase that actually considers the input curves P and 
Q. Our eventual goal is to compute the intersection of reach(P, Q) with the cell boundaries, while 
taking advantage of the data structure from Section[3} For this, we aggregate the cells of FSD(P, Q) 
into (concrete) elementary boxes consisting of r x r cells. There are n 2 /r 2 such boxes. Note that 
we can avoid rounding issues by either duplicating vertices or handling a small part of FSD(P, Q) 
without lookup tables. 

The goal is to determine the signature for each elementary box S. However, at this point this 
is not quite possible yet, because the full signature of S depends on the intersection of reach(i- > , Q) 
with the lower and left boundary of S. Nonetheless, we can find the partial signature, in which the 
positions of so,io (the reach-door) in the (partial) door-orders cr[, aj are still to be determined. 

We aggregate the columns of FSD(P, Q) into vertical strips, each corresponding to a single 
column of elementary boxes (that is, r consecutive columns of cells in FSD(P, Q)). See Figure [7J 

Let A be such a strip. It corresponds to a subcurve P' of P with r edges. The following 
lemma implies that we can build a data structure for A such that, given any segment of Q, we can 
efficiently find its partial door-order within the elementary box in A. 

Lemma 4.1. There exists a constant c, such that the following holds: given a subcurve P' with r 
edges, we can compute in 0{t c ) time a data structure that requires 0{t c ) space and that allows us 
to determine the partial door-order of any line segment on Q in time O(logr). 
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Fig. 



7: The free-space diagram is subdivided into n 2 /r 2 elementary boxes of size r x r. A strip is 
a column of elementary boxes, and it corresponds to a subcurve of P with r edges. 



Proof. Consider the arrangement A of unit circles whose centers are the vertices of P' (see Figure[8]). 
The partial door-order of a line segment s is determined by the intersections of s with the arcs of 
A (and for a circle not intersecting s by whether s lies inside or outside of the circle). Let l s be the 
line spanned by line segment s. Let £ sa be the line parallel to l s that lies above t s and has distance 
1 from l s . Let £ s b be defined similarly, but below £ s . Suppose we wiggle l s (and the corresponding 
f-sai tsb)- Then the order of intersections of l s and the arcs of A changes only when t s moves over 
a vertex of A or if i s leaves or enters a circle. The latter case corresponds to £ sa or £ s b moving over 
the center of a circle. 

Let V be the set of all vertices of A, and let V' be the union of V with the vertices of P'. We 
compute the arrangement B of the lines dual to V . The arrangement B has 0(t 4 ) vertices. We 



build a point location structure for B that is very similar to the structure by Dobkin and Lipton 27 






A 



B 



Fig. 8: By using the arrangement of A, defined by unit circles centered at vertices of P' , we can 
determine the partial door-order of each segment s on Q. This is done by locating the dual 
points of l s , l sa , and l s b in the dual arrangement B. To make this process efficient, the dual 
arrangement is split into so-called slabs. 
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For this, we subdivide B into strips by drawing a vertical line through each vertex of B. Call the 
resulting slab subdivision S. There are 0(r 4 ) slabs, and each slab has 0(r 2 ) cells. 

Now, consider the triple (p(s) of cells in S that contain the dual points £*, £* a , and £* b . The 
cells in (j)(s) lie in the same slab, and they completely determine the combinatorial structure of 
the intersection between £ s and A. Thus, for every possible triple (f)(s) of cells in a strip of S, we 
construct a list £^( s ) that represents the combinatorial structure of t s n A. There are 0(t 4 • r 6 ) = 
0(t 10 ) such lists, each having size O(r). We can compute A/>(s) by traversing the zone of £ s in A. 
Since circles intersect at most twice and also a line intersects any circle at most twice, the zone 



has complexity 0{r2 a ^), where ct(-) denotes the inverse Ackermann function 41 Theorem 5.11]. 
Since 0{r2 a ^) C (3(t 2 ), we can compute all lists in 0(r 12 ) time. 

Given the list Lm 8 \ , the partial door-order of s is determined by the position of the endpoints 
of s in LM s y There are 0(t 2 ) possible ways for this, and we build a table Tm 8 \ that represents 
them. For each entry in T^uy we store a representative for the corresponding partial door-order. 
As described in the previous section, the representative is a positive integer in the word RAM 
model and a pointer to the appropriate record on a pointer machine. 

The total size of the data structure is 0(r 12 ) and it can be constructed in the same time. A 
query works as follows: given s, we can compute £*, £* sa and £* b in constant time. Our data structure 
contains a balanced binary search tree for the point location. The tree is layered: the first layer 
helps us find the vertical slab containing the dual points, the next three layers allow us to determine 
(j>(s), and the last two layers determine the position of the endpoints of s in the list Lm b \. In total, 
there are a constant number of layers of depth O(logr), so we can find a representative for the 
door-order for s in O(logr) time. This bound holds both on the word RAM and on the pointer 
machine. □ 

We can now prove the following lemma, independent of the computational model. 

Lemma 4.2. Given the data structure of Lemma \3.2\ the partial signature for each elementary box 
can be determined in time 0(jit c ~ x + n 2 (log t)/t) for some constant c. 



Proof. By building and using the data structure from Lemma 4.1 we can determine the partial 



door-order for each row in each vertical r-strip in total time proportional to 

n. . , . -i n 2 logr 
-(r c + n log r) = nr c 1 H — . 



T 



We repeat the procedure with the horizontal strips. Now we know for each elementary box 
in FSD(P, Q) the partial door-order for each row and each column. We can then use the data 



structure of Lemma 3.2 to put these together. Since there are n 2 /r 2 boxes, the number of steps is 
0{n 2 /t + n) = 0(n 2 /r). Hence, the partial signature for each elementary box can be computed in 
0{nT c ~ l +n 2 (logr)/r). □ 



5 Solving the Decision Problem 

With the data structures and preprocessing from the previous sections, we have all elements in 
place to determine whether dp(P,Q) < 1. We know for each elementary box its partial signature 
and we have a data structure to derive its full signature (and with it, the combinatorial reachability 
structure) when its reach-doors are known. What remains to be shown is that we can efficiently 



11 



process the free-space diagram to determine whether (n,n) 6 reach(P, Q). This is captured in the 
following lemma. 



Lemma 5.1. If the partial signature for each elementary box is known, we can determine whether 
(n,n) G reach ( P, Q) in time O (n 2 (log t)/t). 

Proof. We go through all of the elementary boxes of FSD(P, Q), processing them one column at a 
time, going from bottom to top in each column. Initially, we know the full signature for the box 
S in the lower left corner of FSD(i- > , Q). We use the signature to determine the intersections of 
reach(P, Q) with the upper and right boundary of S. There is a subtlety here: the signature gives 
us only the combinatorial reachability structure, and we need to map the resulting Si,tj back to 
the corresponding vertices on the curves. On the word RAM, this can be done easily through table 
lookups. On the pointer machine, we use representative records for the Sj, t% elements and use 0(t) 
time before processing the box to store a pointer from each representative record to the appropriate 
vertices on P and Q. 

We proceed similarly for the other boxes. By the choice of the processing order of the elementary 
boxes we always know the incoming reach-doors on the bottom and left boundary when processing 
a box. Given the incoming reach-doors, we can determine the full signature and find the structure 



of the outgoing reach-doors in total time 0(r log r), using Lemma 3.2 Again, we need 0(t) 
additional time on the pointer machine to establish the mapping from the abstract Si, U elements 
to the concrete vertices of P and Q. In total, we spend 0(r log r) time per box. Thus, it takes 
time O (n 2 (log r)/r) to process all boxes, as claimed. □ 

As a result, we obtain the following theorem for a pointer machine (and by extension, for the 
real RAM model). For the word RAM model, we can obtain an even faster algorithm, as described 
in the next section. 

Theorem 5.2. The decision version of the Frechet problem can be solved in time 

0(n 2 (log log n) 3/2 / v^og n) 

on a pointer machine. 

Proof. Set t = A ylog nj log log n for a sufficiently small constant A > 0. The theorem follows by 



applying Lemmas 3.2, 4.2 and 5.1 in sequence. □ 



6 Improved Bound on Word RAM 

We now explain how the running time of our algorithm can be improved if our computational 
model allows for constant time table- lookup. We use the same t as above (up to a constant factor). 
However, we change a number of things. "Signatures" are represented differently and the data 
structure to obtain combinatorial reachability structures is changed accordingly. Furthermore, we 
aggregate elementary boxes into clusters and determine "partial door-orders" for multiple boxes at 
the same time. Finally, we walk the free-space diagram based on the clusters to decide dp(P, Q) < 1. 

Clusters and extended signatures. We introduce a second level of aggregation in the free-space 
diagram (see Figure |9J): a cluster is a collection of r x r elementary boxes, that is, r 2 x r 2 cells 
in FSD(P, Q). Let R be a row of cells in FSD(P, Q) of a certain cluster. As before, the row R 
corresponds to an edge e on Q and a subcurve P' of P with r 2 edges. We associate with R an 
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Fig. 9: A cluster consists of r x r elementary boxes, thus of r 2 x r 2 cells. A row R and its 
corresponding R' for the central elementary box are indicated. 

ordered set Z = (eo, z' , Zi, z[, z%, z' 2 , ■ ■ ■ , Zk, z' k , e\) with 2 • k + 3 elements. Here k is the number 
of intersections of e with the unit circles centered at the r vertices of P 1 (all but the very first). 
Hence, k is bounded by 2r and \Z\ is bounded by 4r + 3. The order of Z indicates the order of 
these intersections with e directed along Q. The elements eo and e\ represent the endpoints of e 
and take a special role. In particular, these elements are used to represent closed doors and snap 
open doors to the edge e. The elements z[ are placeholders for the positions of the endpoints of the 
reach-doors: z' represents a possible reach-door endpoint between eo and z\, the element z[ is an 
endpoint between Z\ and zi, etc. 

Consider a row R' of an elementary box inside the row R of a cluster, corresponding to an edge e 
of Q. The door-index of R' is an ordered set {sq, to, ... , s T , t T ) of size 2r+2. Similar to a door-order, 
elements so and to represent the reach-door at the leftmost boundary of R'; the elements Sj and 
t% (1 < i < t) represent the door at the right boundary of the « th cell in R'. However, instead of 
rearranging the set to indicate relative positions, the elements Sj and ij simply refer to an element 
in Z. If the door is open, then they refer to the intersections with e (possibly snapped to eo or 
ei). If the door is closed, then Si is set to e\ and U is set to eo- Again, the elements so and to are 
special, representing the reach-door and thus refer to one of the elements z[. A partial door-index 
is a door- index without so and to- The advantage of a door- index over a door-order is that the 
reach-door is always at the start. Hence, completing a partial door-index to a full door- index can 
be done in constant time. Since a door-index has size 2r + 2, the number of possible door-indices 
for R! is r°( T ). 

We define the door-indices for the columns analogously, and we concatenate the door-indices 
for the rows and the columns to obtain the indexed signature for an elementary box. Similarly, 
we define the partial indexed signature. The total number of possible indexed signatures remains 

r O(r 2 )_ 

For each possible partial indexed signature S we build a lookup table Xs as follows: the input 
is a word with 4r fields of O(logr) bits each. Each field stores the positions in Z of the endpoints 
of the ingoing reach-doors for the elementary box: 2r fields for the left side, 2r fields for the lower 
side. The output consists of a word that represents the indices for the elements in Z that represent 
the outgoing reach-doors for the upper and right boundary of the box. Thus, the input of Ts is a 
word of 0(t log r) bits, and Is has size t°( t > . Hence, for all partial indexed signatures combined, 
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Fig. 10: (left) Every field represents the partial door-index of a row in an elementary box. (center) 
The fields are grouped into words per row in a cluster, (right) Transposition yields the 
desired organization, where a word represents the partial door-index of the rows in an 
elementary box. 



the size is t°^ t2 ^ = o(n) by our choice of r. 

Preprocessing a given input. During the preprocessing for a given input P, Q, we use superstrips 
consisting of r strips. That is, a superstrip is a column of clusters and consists of r 2 columns of 



the free-space diagram. Lemma 4.1 still holds, albeit with a larger constant c. The data structure 
gets as input a query edge e, and it returns in 0(log7~) time a word that contains r fields. Each 
field represents the partial door-index for e in the corresponding elementary box and thus consists 
of 0(t log r) bits. Hence, the word size is 0(r 2 log r) = O(logn) by our choice of r. Thus, 
the total time for building a data structure for each superstrip and for processing all rows is 
0(n/r 2 (r c + re log r)) = 0(n 2 (log r)/r 2 ). We now have parts of the partial indexed signature for 
each elementary box packed into different words. To obtain the partial indexed signature, we need 
to rearrange the information such that the partial door-indices of the rows in one elementary box 
are in a single word. This corresponds to computing a transposition of a matrix, as is illustrated 
in Figure [lO| For this, we need the following lemma, which can be found — in slightly different 
form — in Thorup |42[ Lemma 9]. 

Lemma 6.1. Let X be a sequence ofr words that contain r fields each, so that X can be interpreted 
as a t x r matrix. Then we can compute in time 0{r log r) on a word RAM a sequence Y of r 
words with r fields each that represents the transposition of X. 

Proof. The algorithm is recursive and solves a more general problem: let X be a sequence of a 
words that represents a sequence M of b different a x a matrices, such that the i th word in X 
contains the fields of the i th row of each matrix in M from left to right. Compute a sequence of 
words Y that represents the sequence M' of the transposed matrices in M. 

The recursion works as follows: if a = 1, there is nothing to be done. Otherwise, we split X 
into the sequence X\ of the first a/2 words and the sequence X<i of the remaining words. X\ and 
X2 now represent a sequence of 26 (a/2) x (a/2) matrices, which we transpose recursively. After 
the recursion, we put the (a/2) x (a/2) submatrices back together in the obvious way. To finish, we 
need to transpose the off-diagonal submatrices. This can be done simultaneously for all matrices 
in time 0(a), by using appropriate bit-operations (or table lookup). 

Hence, the running time obeys a recursion of the form T(a) = 2T(a/2) + 0(a), giving T(a) = 
0(a log a), as desired. □ 

By applying the lemma to the words that represent r consecutive rows in a superstrip, we obtain 
the partial door-indices of the rows for each elementary box. This takes total time 0((n/r 2 ) • (n/r) • 
rlogr) = O (re 2 (log r)/r 2 ). We repeat this procedure for the horizontal superstrips. By using an 
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appropriate lookup table to combine the partial door-indices of the rows and columns, we obtain 
the partial indexed signature for each elementary box in total time 0(n 2 (log t)/t 2 ). 
The actual computation. We traverse the free-space diagram cluster by cluster (recall that a 
cluster consists of r x r elementary boxes). The clusters are processed column by column from 
left to right, and inside each column from bottom to top. Before processing a cluster, we walk 
along the left and lower boundary of the cluster to determine the incoming reach-doors. This is 
done by performing a binary search for each box on the boundary, and determining the appropriate 
elements z[ which correspond to the incoming reach-doors. Using this information, we assemble the 
appropriate words that represent the incoming information for each elementary box. Since there 
are n 2 /r 4 clusters, this step requires time 0((n 2 /t 4 )t 2 log r) = 0(n 2 (log r)/r 2 ). We then process 
the elementary boxes inside the cluster, in a similar fashion. Now, however, we can process each 
elementary box in constant time through a single table lookup, so the total time is 0(n 2 /t 2 ). Hence, 
the total running time of our algorithm is O (n 2 (log r)/r 2 ). By our choice of r = Ay^ogn/log logn 
for a sufficiently small A > 0, we obtain the following theorem. 

Theorem 6.2. The decision version of the Frechet problem can be solved in time 

0(n 2 (log logn) 2 / logn) 

on a word RAM machine. □ 
7 Computing the Frechet Distance 

The optimization version of the Frechet problem, that is, computing the Frechet distance, can be 
done in 0(n 2 logn) time using parametric search with the decision version as a subroutine 16]. We 
showed that the decision problem can be solved in o(n 2 ) time. This however does not directly 
yield a faster algorithm for the optimization problem: if the running time of the decision problem 
is T(n) steps, parametric search results in an 0((T(n) + n 2 )logn) time algorithm j6j. There is 
an alternative randomized algorithm by Raichel and Har-Peled [37], which however also runs in 
0((T(n) + n 2 ) logn) time. We adapt this algorithm to speed up the optimization problem. 

Before we do so, we recall that possible values of the Frechet distance are limited to a certain 
set of critical values [6]: 

1. The distance between a vertex of the one curve and a vertex of the other curve (vertex- 
vertex) , 

2. The distance between a vertex of the one curve and an edge of the other curve (vertex-edge), 

3. For two vertices of one curve and an edge of the other curve, the distance between one of 
the vertices and the intersection of e with the bisector of the two vertices (if this intersection 
exists) (vertex- vertex-edge). 

If we also include vertex-vertex-edge tuples with no intersection, we can sample a critical value 
uniformly at random in constant time. The algorithm now works as follows (see Har-Peled and 
Raichel [37] for more details): it first samples K = <in 2 critical values uniformly at random. Next 
the algorithm finds the interval [a, b] with a and b being two critical values in the sample such that 
the Frechet distance lies in [a, b], and no other critical value of the sample lies in [a, b]. In the 
original algorithm this is done by sorting the critical values and performing a binary search using 
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the decision version. By using median- finding instead, this step can be done in 0{K + T(n) log K) 
time. We note that, alternatively, the running time of this step could be reduced by picking a 
smaller K. However, this does not improve the final bound, since it is dominated by a 0(n 2 2 a ^) 
component. 

This atomic interval [a, b] of the sampled critical values with high probability contains only a 
small number of the remaining critical values. More specifically, for K = An 2 the probability that 
the interval contains more than 2cn\an critical values at most l/n c |37| Lemma 6.2]. 

The remainder of the algorithm first determines the K' critical values in the interval [a, b], then 
again sorts them and performs a binary search. Excluding the time to determine the critical values, 
this takes 0(K' + T(n) \ogK') time with mean-finding. Thus the crucial part is to determine the 
K' critical values fast. 

In 0(n 2 ) time we can check for any vertex- vertex and vertex-edge pair whether the corre- 
sponding critical value lies in [a, b]. It remains to determine the critical values corresponding 
to vertex-vertex-edge tuples. These critical values are determined by a variant of the standard 
sweepline algorithm. For this take an edge e of P and the vertices of Q. The sweep starts with 
circles of radius a around the vertices of Q and increases the radii until they reach b. During this 
sweep the algorithm maintains the order in which the circle arcs intersect e. A critical value of 
the vertex-vertex-edge type corresponds to the event that two different circles intersect e in the 
same point. Next to these events the sweepline algorithm requires the following events: a circle 
intersects e for the first time, or a circle intersects one of the vertices of e. Both of these event 
types correspond to critical values involving e or a vertex of e. Thus if we perform such a sweep 
for all edges of P (and similarly for the edges of Q), the total number of events is O(K'), thus the 
overall running time of all sweeps ignoring the time for initialization is 0(K' log n). 

It remains to show that we can compute the initial order in which the circle arcs intersect e 
fast. First compute the arrangement A of circles with radius a around the vertices of Q. This 



can be done in 0(n 2 ) time 23 . We need to determine in which order the arcs of the circles 
intersect e. We can determine the order of intersections by traversing in A the zone of the line £ 
spanned by e. The time for the traversal can be bounded by the complexity of the zone. Using 
that the circles pairwise intersect at most twice and the line intersects each circle also only twice, 



the complexity of the zone can be bounded by 0{n2 a ^) (4ll Theorem 5.11]. Summing over all 
edges e this adds a total of 0(n 2 2 a ^) to the running time. Thus the overall running time is 
0{T(n) log(re) + n 2 2 a ^ + K' log n) . The case that K' > 8n In n happens with probability less than 
1/n 4 , and also in this case K' is still in 0(n 3 ). Thus, this case adds o(l) to the expected running 
time. The case K' < 8nlnn adds 0(nlog 2 n) to the expected running time. As a consequence we 
obtain the following lemma. 

Lemma 7.1. The Frechet distance of two polygonal curves with n vertices each can be computed by 
a randomized algorithm in 0(n 2 2 a ^ +T{n) log n) expected time, where T(n) is the running time 
for the decision problem. 

We plug in our new bounds on T(n). 

Theorem 7.2. The Frechet distance of two polygonal curves with n vertices each can be computed 
by a randomized algorithm in time 0(n 2 \/log n(log log n) 3 / 2 ) on a pointer machine and in time 
0(n 2 (loglogn) 2 ) on a word RAM. 
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8 Decision Trees 



Our results also have implications for the decision-tree complexity of the Frechet problem. Since in 
that model we account only for comparisons between the input elements, the preprocessing comes 
for free, and hence the size of the elementary boxes can be increased. 

Before we consider the continuous Frechet problem, we first note that the techniques of Agar- 
wal et al. can be used to obtain a similar result for the discrete Frechet problem: suppose we have 
two sequences P = (pi,P2, ■ ■ ■ ,Pn) and Q = (qi, q2, . ■ ■ , q n ). For 5 > 0, we define a directed graph 
G$ with vertex set P x Q. In G$, there is an edge between two vertices (pi,qj), (pi,Qj+i) if and 
only if both d(pi, qj) < 5 and d(pi, qj+i) < S. The condition is similar for an edge between vertices 
(pi,qj) and (pi + \,qj), and vertices (pi,qj) and (pi+i, Qj+i)- There are no further edges in G$- Now 
the problem is to find the smallest 5 for which G$ has a path from (pi, qi) to (p n , q n ). 

Using the techniques from Agarwal et al. we can obtain the following theorem^] 

Theorem 8.1. There is an algebraic computation tree for the discrete Frechet problem of depth 
0(n 3 / 2 ). 

Proof. We first discuss the decision problem, where we are given P, Q and a 5, and we need to 
decide whether we can reach (p n ,q n ) from (pi,gi). For the discrete case, the analogue of the 
reachable free-space is just an n x n boolean matrix M, where the bit in the i th row and the j th 
column indicates whether the pair (pi,qj) can be reached from (pi,qi) in G&. 

We set t = \fn and subdivide the columns of M into strips of width r, as above. Each strip 
corresponds to a subsequence Q' of r points qj , and we compute the arrangement of unit disks with 
centers in Q' . This takes 0(r 2 ) time. Then we locate each point of P in this arrangement, in total 
time 0(n logn). We do this for every strip, resulting in a total running time of 0(n 3//2 ), by our 
choice of r. As observed by Agarwal et al., the information we gain in this way suffices to complete 
M without further comparisons. As shown by Agarwal et al., one can then solve the optimization 
problem at the cost of another 0(logn)-factor, which is absorbed into the O-notation. □ 

Given our results from above, we can prove an analogous statement for the continuous Frechet 
distance. 

Theorem 8.2. There exists an algebraic decision tree for the Frechet problem (decision version) 
of depth 0(n 2 ~ £ ), for a fixed constant e > 0. 

Proof. We reconsider the steps of our algorithm. The only phases that actually involve the input are 
the second preprocessing phase and the traversal of the elementary boxes. The reason of our choice 
for t was to keep the time for the first preprocessing phase polynomial. This is no longer a problem. 



By Lemmas 4.2 and 5.1 the cost for t he r emaining algorithm is bounded by 0(nr° +n (log t) /t) 



where c is the constant from Lemma 4.1 Choosing r = n 1//c , we get a decision tree of depth 



r i ir log /> / t , , , 

n ■ n^r -| -f— = 0(n 2 ~ [1/c > log n) = 0{n 



2-e\ 

— > p \ - I ill ; i ) I fill 



for, say, e = l/2c. □ 
Assuming linear time reductions, this has the following consequence for Alt's conjecture: 
8 The notation O(-) (pronounced "soft-Oh") stands for O(-) up to poly-logarithmic factors. 
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Corollary 8.3. // the decision version of the Frechet problem is 3S\JM-hard, then 3SUM has an 
algebraic decision tree of depth 0(n 2 ~ 6 ). 

We leave it to the reader to judge the implications on the status of the conjecture. 
9 Conclusion 

In this paper we have broken the long-standing quadratic upper bound for the decision version of 
the Frechet problem. Moreover, we have shown that this problem has an algebraic decision tree 
of depth 0(n 2_£ ), for some e > and where n is the number of vertices of the polygonal curves. 
This strongly indicates that the problem is not 3SUM hard after all. We have shown how our 
faster algorithm for the decision version can be used for a faster algorithm to compute the Frechet 
distance. If we allow constant-time table-lookup, we obtain a running time in close reach of 0(n 2 ). 

This leaves us with intriguing open research questions. Can we reduce the time needed for the 
decision version to 0{n 2 ~ e ), the bound we obtain from the algebraic decision tree? Can we devise 
a quadratic or even subquadratic algorithm for the optimization version? Can we devise such an 
algorithm on the word RAM, that is, with constant-time table-lookup? Or, on the other hand, 
can we establish a connection between the Frechet distance and other problems which exhibit 
a discrepancy between the decision tree and the uniform complexity, such as, e.g., Min-Plus- 
Convolution? 
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A Computational Models 

Real RAM.. The standard machine model in computational geometry is the real RAM. Here, 
data is represented as an infinite sequence of storage cells. These cells can be of two different types: 
they can store real numbers or integers. The model supports standard operations on these numbers 
in constant time, including addition, multiplication, and elementary functions like square-root, sine 
or cosine. Furthermore, the integers can be used as indices to memory locations. Integers can be 
converted to real numbers in constant time, but we need to be careful about the reverse direction. 
The floor function can be used to truncate a real number to an integer, but if we were allowed 
to use it arbitrarily, the real RAM could solve PSPACE-complete problems in polynomial time. 
Therefore, we usually have only a restricted floor function at our disposal. 

Word RAM.. The word RAM is essentially a real RAM without support for real numbers. 
However, on a real RAM, the integers are usually treated as atomic, whereas the word RAM 
allows for powerful bit-manipulation tricks. More precisely, the word RAM represents the data as 
a sequence of u>-bit words, where w = r2(logn). Data can be accessed arbitrarily, and standard 
operations, such as Boolean operations (and, xor, shl, . . .), addition, or multiplication take constant 
time. There are many variants of the word RAM, depending on precisely which instructions are 
supported in constant time. The general consensus seems to be that any function in AC is 
acceptable^] However, it is always preferable to rely on a set of operations as small, and as non- 
exotic, as possible. Note that multiplication is not in AC , but nevertheless is usually included in 
the word RAM instruction set. 

9 AC is the class of all functions /: {0, 1}* — > {0, 1}* that can be computed by a family of circuits (C„)„gM with 
the following properties: (i) each C„ has n inputs; (ii) there exist constants a, b, such that C n has at most an b gates, 
for n6N; (iii) there is a constant d such that for all n the length of the longest path from an input to an output in 
Cn is at most d (ie, the circuit family has bounded depth) ; (iv) each gate has an arbitrary number of incoming edges 
(ie, the fan-in is unbounded). 
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Pointer machine.. The pointer machine model disallows the use of constant time table lookup, 
and is therefore a restriction of the (real) RAM model. The data structure is modeled as a directed 
graph G with bounded out-degree. Each node in G represents a record, with a bounded number 
of pointers to other records and a bounded number of (real or integer) data items. The algorithm 
can access data only by following pointers from the inputs (and a bounded number of global entry 
records); random access is not possible. The data can be manipulated through the usual real RAM 
operations, but without support for the floor function, for reasons mentioned above. 

Algebraic computation tree.. Algebraic computation trees (ACTs) [9j are the computational 
geometry analogue of binary decision trees, and like these they are mainly used for proving lower 
bounds. Let x\,... ,x n £ K be the inputs. An ACT is a binary tree with two different kinds of 
nodes: computation nodes and branch nodes. A computation node v has one child and is labeled 
with an expression of the type y v = y u @ y w , where © G {+,—,*,/, \f-} is a operation and y u , y w is 
either an input variable x\, . . . ,x n or corresponds to a computation node that is an ancestor of v. 
A branch node has degree 2 and is labeled by y u = or y u > 0, where again y u is either an input 
or a variable corresponding to an ancestor. A family of algebraic computation trees (T n ) ng N solves 
a computational problem (like Delaunay triangulation or convex hulls computation), if for each 
r£N, the tree T n accepts inputs of size n, and if for any such input x±, . . . ,x n the corresponding 
path in T n (where the children of the branch nodes are determined according the conditions they 
represent) constitutes a computation which represents the answer in the variables y v encountered 
during the path. 
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